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Abstract — In measurements of in-room radio channel re- 
sponses an avalanche effect can be observed: earliest signal 
components, which appear well separated in delay, are followed 
by an avalanche of components arriving with increasing rate of 
occurrence, gradually merging into a diffuse tail with exponen- 
tially decaying power. We model the channel as a propagation 
graph in which vertices represent transmitters, receivers, and 
scatterers, while edges represent propagation conditions between 
vertices. The recursive structure of the graph accounts for the 
exponential power decay and the avalanche effect. We derive 
a closed form expression for the graph's transfer matrix. This 
' expression is valid for any number of interactions and is 
straightforward to use in numerical simulations. We discuss an 
example where time dispersion occurs only due to propagation in 
between vertices. Numerical experiments reveal that the graph's 
recursive structure yields both an exponential power decay and 
an avalanche effect. 

I I. Introduction 

Engineering of modem indoor radio systems for communi- 
cations and geolocation relies heavily on models for the time 
dispersion of the wideband and ultrawideband radio channels 
. Q-lUl. From measurement data, as exemplified in Fig. [T] 
■ it appears that the spatial average of the channel impulse 
, response's squared magnitude (referred to as the delay-power 
' spectrum) observed in in-room scenarios exhibits an avalanche 
effect: The earliest signal components, which appear well 
separated in time, are followed by an avalanche of components 
, arriving with increasing rate of occurrence, gradually merging 
' into a diffuse tail with exponentially decaying power The 
diffuse tail is often referred to as a "dense multipath compo- 
nent" and is commonly attributed to diffuse scattering from 
rough surfaces and objects which are small compared to the 
wavelength [5J. A similar avalanche effect is well-known in 
, room acoustics |l6l where it is attributed to recursive scattering 
of sound waves. Indoor radio propagation environments are 
particularly exposed to recursive scattering as electromagnetic 
waves may be reflected back and forth in between walls, floor, 
and ceiling. Thus, the avalanche effect and the exponential 
power decay may occur due to recursive scattering rather than 
diffuse scattering. 

Recursive scattering phenomena have been previously taken 
into account in a number of radio channel models. The works 
iTTl- lfTOl use the analogy to acoustical reverberation theory to 
predict the exponential decay. As a matter of fact, there exists 
a well-developed theory of electromagnetic fields in cavities 
IfTTI . ||T21 . but in this context too the avalanche effect has 
received little attention. Recursive scattering between particles 
in a homogeneous medium is a well-known phenomenon 
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Fig. 1. Spatially averaged delay-power spectrum measured within an office 
of 5 x5x2. 6 m"^ obtained as the rms value of impulse response for 30x30 
receiver positions on a square horizontal grid with steps of 1 cm. Abscissa 
includes cable delays; reference to signal bandwidth (10 GHz) shifts ordinate 
by -100 dB. Reprinted from |4) wifli permission (© 2002 IEEE). 



Studied by Foldy fT3\ and Lax (14], (15] . The solution, given 
by the so-called Foldy-Lax equation |16|, has been applied 
in the context of time-reversal imaging by Shi and Nehorai 
p7l. The solution is, however, intractable for heterogeneous 
indoor environments. In [18] the radio propagation mechanism 
is modeled as a "stream of photons" performing continuous 
random walks in a homogeneously cluttered environment. The 
model predicts a delay-power spectrum consisting of a single 
directly propagating "coherent component" followed by an 
incoherent tail. Time-dependent radiosity |19|-|22| accounting 
for delay dispersion has been recently applied to design a 
model for the received instantaneous power [23]. Thereby, 
the exponential power decay and the avalanche effect can be 
predicted. 

Simulation studies of communication and localization sys- 
tems commonly rely on synthetic realizations of the channel 
impulse response. A multitude of impulse response models 
exist Q, 12], ll24l . but only few account for the avalanche 
effect. The models Il25l - ll27l treat early components via a ge- 
ometric model whereas the diffuse tail is generated via another 
stochastic process; the connection between the propagation 
environment and the diffuse tail is, however, not considered. 
Ray tracing methods may also be used to do site-specific 
simulations ll28l . However, ray tracing methods commonly 
predict the signal energy to be concentrated into the single 
dominant components whereas the diffuse tail is not well 
represented. 

In this contribution, we model the channel impulse response 
for the purpose of studying the avalanche effect and, in 
particular, its relation to recursive scattering. The objective 
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is a model which allows for simulation of both realizations of 
the channel response and average entities such as the delay- 
power spectrum. Expanding on the previously presented work 
||29l . OOll . we propose a unified approach for modeling of 
the transition from single dominant components to the diffuse 
tail. The propagation environment is represented in terms of 
a propagation graph, with vertices representing transmitters, 
receivers, and scatterers, while edges represent propagation 
conditions between vertices. The propagation graph accounts 
inherently for recursive scattering and thus signal components 
arriving at the receiver may have undergone any number 
of scatterer interactions. This modelling approach allows for 
expressing the channel transfer function in closed form for 
unlimited number of interactions. The formalism enables a 
virtual "dissection" of the impulse response in so-called partial 
responses to inspect how recursive scattering leads to a gradual 
buildup of the diffuse tail. 

Propagation graphs may be defined according to a spe- 
cific scenario for site-specific prediction, or be generated 
as outcomes of a stochastic process for use in e.g. Monte 
Carlo simulations. In the present contribution we consider 
an example of a stochastically generated propagation graph 
suitable for Monte Carlo simulations. In the example model, 
scatterer interactions are assumed to cause no time dispersion 
and thus delay dispersion occurs only due to propagation 
in between vertices. The simulations reveal that the graph's 
recursive structure yields both an exponential power decay and 
an avalanche effect in the generated impulse responses. 

II. Representing Radio Channels as Graphs 

In a typical propagation scenario, the electromagnetic signal 
emitted by a transmitter propagates through the environment 
while interacting with a number of objects called scatter- 
ers. The receiver, which is usually placed away from the 
transmitter, senses the electromagnetic signal. If a line of 
sight exists between the transmitter and the receiver, direct 
propagation occurs. Also, the signal may arrive at the receiver 
indirectly via one or more scatterers. In the following we 
represent propagation mechanisms using graphs allowing for 
both recursive and non-recursive scattering. We first introduce 
the necessary definitions of directed graphs. 

A. Directed Graphs 

Following II3TI . we define a directed graph Q as a pair 
(V, £) of disjoint sets of vertices and edges. Edge e <E £ with 
initial vertex denoted by init(e) and terminal vertex denoted by 
term(e) is said to be outgoing from vertex init(e) and ingoing 
to term(e). We consider graphs without parallel edges. Thus, 
there exists no pair of edges e and e' such that init(e) = 
init(e') and term(e) = term(e'). In this case an edge e can 
be identified by the vertex pair (init(e), term(e)) £ and 
with a slight abuse of notation we write e = (init(e), term(e)) 
and f C V X V. 

B. Propagation Graphs 

We define a propagation graph as a directed graph where 
vertices represent transmitters, receivers and scatterers. Edges 




£d = {(Txl, Rxl), (Tx2, Rxl), (Tx3, Rx;l)} 

£, = {(Tx2,S3),(Tx3,S3),(Tx4,S6),(Tx4,Sl)} 

= {(SI, Rx3), (S3, Rx3), (S6, Rx2)} 
£a = {(SI, S2), (S2, SI), (S3, SI), (S2, S4), (S4, S3), (S4. S5)} 



Fig. 2. A propagation graph witli four transmit vertices, three receive vertices 
and six scatterer vertices. 

represent the propagation conditions between the vertices. 
Thus, the vertex set V of a propagation graph is a union 
of three disjoint sets: V = Vt U Vr U Vs, the set of 
transmitters Vt — {Txl, . . . , TxiVt}, the set of receivers 
Vr = {Rxl, . . . , RxA'i}, and the set of scatterers Vs = 
{SI, ... , SiVg}. The transmit vertices are considered as sources 
and have only outgoing edges. Likewise, the receivers are 
considered as sinks with only incoming edges. The edge set 
£ can thus be partitioned into four disjunct sets as £ = 
£dU £tU £i- U £s, with direct edges in = ^ n (Vt x Vr), 
transmitter-scatterer edges in ft = ^ H (Vt x Vs), scatterer- 
receiver edges in £i- = £0 (Vs x Vr), and inter-scatterer edges 
in £s — £ n (Vs X Vs). Fig. |2] shows an example propagation 
graph. 

The signals propagate in the graph in the following way. 
Each transmitter emits a signal that propagates via its outgoing 
edges. The signal observed by a receiver vertex is the sum 
of the signals arriving via the ingoing edges. A scatterer 
sums up the signals on its ingoing edges and re-emits the 
sum-signal on the outgoing edges. As a signal propagates 
along an edge, or interacts with a scatterer, it undergoes 
delay and dispersion in time. The specific delay dispersion 
endured by a signal depends on the particular propagation 
mechanism along its edges. Assuming these mechanisms to 
be linear and time-invariant, this effect can be represented 
as a convolution with an impulse response or, in the Fourier 
domain, as a multiplication with a transfer function. Hence, 
the signal arriving to vertex u„/ via edge e = (w„,u„') reads 
Ae(/)C„(/), where C„(/) is the signal emitted by vertex u„ 
and Ae{f) denotes the transfer function associated to edge 
e. In other words, the transfer function Ae{f) describes the 
interaction at the initial vertex u„ and the propagation from 

Vn to Vn'- 

C. Weighted Adjacency Matrix of a Propagation Graph 

Propagation along the edges is described via a transfer 
matrix A(/) which can be viewed as an adjacency matrix 
with each entry coinciding with the edge transfer function of 
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Fig. 3. Vector signal flow grapli representation of a propagation graph. 
Vertices represent vertex sets of the propagation graph with associated vector 
signals. Signal transmission between the sets are represented by the edges and 
associated transfer matrices. 



the corresponding edge. Thus, the weighted adjacency matrix 
A(/) G C^N,+N^+N^)x(N,+N^+N^) of jhe propagation graph 

G is defined as 



[A(/)]. 



^(«„, «„/)(/) if («n,w„') e £, 







otherwise. 



(1) 



i.e., entry n,n' of A(/) is the transfer function from vertex 
Vn to vertex u„/ of Q. Selecting the indexing of the vertices 
according to 



v„ G 



Vt, 

Vr, 



n — 1, . 

n = Nt 
n^Nt 



l,...,iVt 
N, + 1, . . 



■iVr 

,Nt + N, + Ns, 



(2) 



the weighted adjacency matrix takes the form 

A(/) = 





D(./) R(/) 
T(/) B(/) 



(3) 



where denotes the all-zero matrix of the appropriate dimen- 
sion and the transfer matrices 



D(/)G 

T(/)e 
B(/)e 



connecting Vt to Vr, (4) 

connecting Vs to Vr, (5) 

connecting Vt to Vs, and (6) 

interconnecting Vg. (7) 



The special structure of A(/) reflects the structure of the 
propagation graph. The first iVt rows are zero because we 
do not accept incoming edges into the transmitters. Likewise 
columns A't + 1 , . . . , A't + A'r are all zero since the receiver 
vertices have no outgoing edges. 

The input signal vector X(/) is defined as 



X{f) = [X,{f),...,XN,{f)f, 



(8) 



where Xm{f) is the spectrum of the signal emitted by trans- 
mitter Txrn, and [ J^ denotes the transposition operator The 
output signal vector Y(/) is defined as 



where Y,n{f) is the Fourier transform of the signal observed 
by receiver Rxm. Similar, to X(/) and Y(/) we let Z(/) 
denote the output signal vector of the scatterers: 



z{f) = [z,{f),...,z^M)]' 



(10) 



with the nth entry denoting the Fourier transform of the 
signal observed at scatterer vertex Sn. By the definition of the 
propagation graph, there are no other signal sources than the 
vertices in Vt. Assuming linear and time-invariant propagation 
mechanisms, the input-output relation in the Fourier domain 
reads 

Y(/) =H(/)X(/), (11) 

where H(/) is the N^. x A't transfer matrix of the propagation 
graph. 

The structure of the propagation graph unfolds in the vector 
signal flow graph depicted in Fig.|3] The vertices of the vector 
signal flow graph represent the three sets Vt, Vr, and Vs 
with the associated signals X(/),Y(/), and Z(/). The edge 
transfer matrices of the vector signal flow graph are the sub- 
matrices of A(/) defined in (|4li-(l7]i. 

III. Transfer Matrix of a Propagation Graph 

In the following we derive the transfer matrix of a prop- 
agation graph. In Subsection IIII-AI we first discuss how the 
response of a graph is composed of signal contributions prop- 
agating via different propagation paths. This representation 
is, albeit intuitive, impractical for computation of the transfer 
matrix of graphs with cycles. Thus in Subsections IIII-BI and 
IIII-CI we give the transfer matrix and partial transfer matrices 
of a general propagation graph in closed form. Subsection 
IIII-DI treats the graphical interpretation of reciprocal channels. 
The section concludes with a discussion of related results in 
the literature. 

A. Propagation Paths and Walks 

The concept of a propagation path is a corner stone in mod- 
eling multipath propagation. In the literature, this concept is 
most often defined in terms of the resulting signal components 
arriving at the receiver A shortcoming of this definition is 
that it is often hard to relate to the propagation environment. 
The graph terminology offers a convenient alternative. A walk 
£ (of length K) in a graph Q ~ (V, £) is defined as a 



sequence £ ~ {v^-^\v^^\. 



.,v 



) of vertices in V such 



Y(/) = [Fi (/),..., IV. (/)] 



(9) 



that e£,k = l,...K.We say that ^ is a walk 

from w'^^ to where u^^^ is the start vertex and 

is the end vertex; v'^^\ . . . , u*^^' are called the inner vertices 
of £. We define a propagation path as a walk in a propagation 
graph from a transmitter to a receiver Consequently, all (if 
any) inner vertices of a propagation path are scatters. A signal 
that propagates along propagation path £ traverses K + 1 (not 
necessarily different) edges and undergoes K interactions. We 
refer to such a propagation path as a /iT -bounce path. The 
zero-bounce path £ — {v,v') is called the direct path, from 
transmitter v to receiver v'. As an example, referring to the 
graph depicted in Fig. |2] it is straightforward to verify that 
£i = (Txl,Rxl) is a direct path, ^ = (Tx4, S6, Rx2) is 
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a single-bounce path, and £3 — (Tx4, SI, S2, SI, Rx3) is a 
3-bounce path. 

We denote by L^v' the set of propagation paths in Q 
from transmitter v to receiver v' . The signal received at v' 
originating from transmitter v is the superposition of signal 
components each propagating via a propagation path in L^^' ■ 
Correspondingly, entry of H(/) reads 



where Hi{f) is the transfer function of propagation path i. 

The number of terms in (fTZt equals the cardinality of Cw' , 
which may, depending on the structure of the graph, be finite 
or infinite. As an example, the number of propagation paths 
is infinite if v and v' are connected via a directed cycle in Q, 
i.e. if a propagation path contains a walk from an inner vertex 
back to itself. The graph in Fig. |2] contains two directed cycles 
which are connected to both transmitters and receivers. 

In the case of an infinite number of propagation paths, 
computing Hyy'{f) directly from (fTZt is infeasible. This 
problem is commonly circumvented by truncating the sum in 
(fT2l i to approximate H^yi (/) as a finite sum. This approach, 
however, calls for a method for determining how many terms 
of the sum should be included in order to achieve reasonable 
approximation. 

In the frequently used "ii'-bounce channel models", propa- 
gation paths with more than K interactions are ignored. This 
approach is motivated by the rationale that at each interaction, 
the signal is attenuated, and thus terms in (fT2l l resulting from 
propagation paths with a large number of bounces are weak 
and can be left out as they do not affect the sum much. 
This reasoning, however, holds true only if the sum of the 
components with more than K interactions is insignificant, 
which may or may not be the case. From this consideration, 
it is clear that the truncation criterion is non-trivial as it 
essentially necessitates computation of the whole sum before 
deciding whether a term can be ignored or not. 

B. Transfer Matrix for Recursive and Non-Recursive Propa- 
gation Graphs 

As an alternative to the approximation methods applied to 
the sum (fT2l i we now give an exact closed-form expression 
for the transfer matrix H(/). Provided that the spectral radius 
of B(/) is less than unity, the expression holds true for any 
number of terms in the sum (fT2l i and thus holds regardless 
whether the number of propagation paths is finite or infinite. 

Theorem 1: If the spectral radius of B(/) is less than unity, 
then the transfer matrix of a propagation graph reads 



H(/) - D(/) + R(/)[I - B(/)]-iT(/). 



(13) 



According to Theorem [T] the transfer matrix H(/) consists 
of the two following terms: D(/) representing direct propa- 
gation between the transmitters and receivers and R(/)[I — 
B(/)]~^T(/) describing indirect propagation. The condition 
that the spectral radius of B(/) be less than unity implies that 
for any vector norm || • ||, ||Z(/)|| > ||B(/)Z(/)|| for non- 
zero ||Z(/)||, cf. [32|. For the Euclidean norm in particular 



this condition implies the sensible physical requirement that 
the signal power strictly decreases for each interaction. 

Proof: Let Hfc(/) denote the transfer matrix for all k- 
bounce propagation paths, then H(/) can be decomposed as 



h(/) = Eh,(/), 



A;=0 



(12) where 







lR(/)B'=-i(/)T(/), fc>0 



Insertion of (flST l into (fl4l i yields 

H(/) = D(/) + R(/) ( £ B'~\f)] T(/). 



(14) 



(15) 



(16) 



The infinite sum in (fTSI l is a Neumann series converging to 
[I — B(/)]^^ if the spectral radius of B(/)) is less than unity. 
Inserting this in ( fTSI l completes the proof. ■ 

The decomposition introduced in (fT4l i makes the effect of 
the recursive scattering directly visible. The received signal 
vector is a sum of infinitely many components resulting from 
any number of interactions. The structure of the propagation 
mechanism is further exposed by (fT6t where the emitted vector 
signal is re-scattered successively in the propagation environ- 
ment leading to the observed Neumann series. This allows 
for modeling of channels with infinite impulse responses by 
expression ( fTST l. 

It is possible to arrive at (fTsT i in an alternative, but less 
explicit, manner: 

Proof: It is readily observed from the vector signal flow 
graph in Fig. |3]that Z(/) can be expressed as 



Z(/) = T(/)X(/)+B(/)Z(/). 



(17) 



Since the spectral radius of B(/) is less than unity we obtain 
for Z(/) the solution 



Z(/) = [I-B(/)]-iT(/)X(/). 



(18) 



Furthermore, according to Fig. |3]the received signal is of the 
form 



Y(/)-D(/)X(/)+R(/)Z(/). 



(19) 



Insertion of (fTsT l in this expression yields (fT?t . ■ 
We remark that the above two proofs allow for propagation 
paths with any number of bounces. This is highly preferable, as 
the derived expression (fT3T l is not impaired by approximation 
errors due to the truncation of the series into a finite number 
of terms as it occurs when using i^-bounce models. 

A significant virtue of the expression (fljT l is that prop- 
agation effects related to the transmitters and receivers are 
accounted for in the matrices D(/), T(/) and R(./), but do 
not affect B(/). Consequently, the matrix [I — B(/)]^^ only 
needs to be computed once even though the configuration of 
transmitters and receivers changes. This is especially advanta- 
geous for simulation studies of e.g. spatial correlation as this 
leads to a significant reduction in computational complexity. 



5 



C. Partial Transfer Matrices 

The closed form expression (fTjt for the transfer matrix of 
a propagation graph accounts for propagation via an arbitrary 
number of scatterer interactions. For some appHcations it is, 
however, relevant to study only some part of the impulse 
response associated with a particular number of interactions. 
One case is where a propagation graph is used to generate 
only a part of the response and other techniques are used for 
the remaining parts. Another case is when one must assess 
the approximation error when the infinite series is truncated. 
In the following we derive a few useful expressions for such 
partial transfer matrices. 

We define the K : L partial transfer matrix as 



H 



K:L 



if) 



L 

E 

k=K 



Hfc(/), < X < L, 



(20) 



i.e., we include only contributions from propagation paths with 
at least K, but no more than L bounces. It is straightforward 
to evaluate ^ for X = 0, and L = 0, 1, 2: 



Ho:o(/)=D(/) 
Ho:l(/)=D(/) 
Ho:2(/)=D(/) 



R.(/)T(/) 

R(/)T(/) 



R(/)B(/)T(/). 



(21) 
(22) 
(23) 



This expansion of the truncated series is quite intuitive but 
the obtained expressions are increasingly complex for large 
L. Theorem |2] gives a closed form expression of the partial 
transfer function Ha',l(/) for arbitrary K and L: 

Theorem 2: The partial response 'H.K-.Lif ) is given by 



H 



K:L 



if) 



D(/) + R(/)[I - B^(/)][I - B(/)]-iT(/), 

K = 0,L>0 

R(/)[B^'-i(/) - B^(/)][I - B(/)]-iT(/), 

0<K <L, 



provided that the spectral radius of B(/) is less than unity. 
Proof: The partial transfer function for < K < L reads 



k=K k'=L + l 

= Hi<-:oo(/) - Hl+1:oc(/)- 



(24) 



For K = we have Ho:oo(/) — H(/) by definition; for 
K > 1 we have 



Hx:oo(/)=R(/) J2 B'=(/)T(/) 

oo 

= R(/)B^-i(/)^B'=(/)T(/) 



fe=0 



- R(/)B^-i(/)[I - B(/)]-iT(/). (25) 

Inserting (l25T l into (l24l l completes the proof. ■ 
Theorem |2] enables closed-form computation of H.K:L{f) 
for any K > L. We have already listed a few partial transfer 
matrices in (ISTT l. (l22l i. and (|23] |. By definition the partial 
response Ha':a:(/) equals Ha-(/) for which an expression 
is provided in ( fTSl l. The transfer function of the iC-bounce 



approximation is equal to Ho:A-(/)- Another special case 
worth mentioning is Hk+i-.ooH) = H(/)— Ho:A-(/) available 
from (l25T l. which gives the error due to the iiT -bounce approx- 
imation. Thus the validity of the X-bounce approximation 
can be assessed by evaluating some appropriate norm of 

Ha-+1:oo(/)- 



D. Reciprocity and Reverse Graphs 

In most cases, the radio channel is considered reciprocal. As 
we shall see shortly, the graph terminology accommodates an 
interesting interpretation of the concept of reciprocity. For any 
propagation graph we can define the reverse graph in which 
the roles of transmitter and receiver vertices are swapped. 
The principle of reciprocity states that the transfer matrix 
of the reverse channel is equal to the transposed transfer 
matrix of the forward channel, i.e., a forward channel with 
transfer matrix H(/) has a reverse channel with transfer matrix 
H(/) — H^(/). In the sequel we mark all entities related to 
the reverse channel with a tilde. 

We seek the relation between the forward graph Q = (V, £) 
and its reverse Q — (V, £) under the assumption of reciprocity. 
More specifically, we are interested in the relation between 
the weighted adjacency matrix A(/) of Q and the weighted 
adjacency matrix A(/) of Q. We shall prove the following 
theorem: 

Theorem 3: For a propagation graph Q = (V, £) with 
weighted adjacency matrix A(/) and transfer matrix H(/) 
there exists a reverse graph Q = (V, £) such that V = V, £ = 
{(«,«') : e £}, A(/) = AT(/), and H(/) = HT(/). 

Proof: We prove the Theorem by constructing a suit- 
able propagation graph such that the reciprocity condition is 
fulfilled. We first note that the set of transmitters, receivers, 
and scatterers is maintained for the reverse channel, thus the 
vertex set of ^ is V with Vt = K, K = Vt, and Vg = Vg. It 
is immediately clear from the structure of Q that the reverse 
graph Q has no ingoing edges to H and no outgoing edges 
from Vi and Q is thus a propagation graph. Assuming the 
vertex indexing as in (|2]l, the weighted adjacency matrix of Q 
is of the form 



A(/) = 



D(/) T(/)- 


R(/) B(/). 



(26) 



where the transfer matrices D(/), R(/), T(/), and B(/) are 
defined according to Fig. |4] Equating A(/) and A(/) we 
obtain the identities D(/) = DT(/), B(/) = B^{f),'T{f) = 
B7{f), and R(/) = T^{f). The relation between £ and £ 
now follows from the definition of the weighted adjacency 
matrix. The input-output relation of the reverse channel reads 
Y(/) = H(/)X(/) where X(/) is the signal vector emitted 
by the vertices in Vt and Y(/) is the signal vector received 
by the vertices in Vt. Considering Fig. |4] and arguing as in 
Section IIII-BI yields for the reverse channel 



H(/) = D(/)+R(/) I-B(/) t(/) 



(27) 
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Fig. 4. Vector signal flow graph representation of a reverse propagation graph. 
Compared to the forward graph depicted in Fig. |3] all edges are reversed. 

Inserting for D(/), R(/), f (/), and B(/) leads to 

H(/) = DT(/) + tT(/) [I - bT(/)] ' rT(/) 

= HT(/), (28) 

where the last equality results from (ITjt . ■ 
In words Theorem |3] says that for a propagation graph there 
exists a reverse graph which fulfills the reciprocity condition. 
Furthermore, an example of a propagation graph fulfilling 
the reciprocity condition is the reverse graph Q obtained by 
reversing all edges of Q while maintaining the edge transfer 
functions. 

E. Related Recursive Scattering Models 

We provide a few examples of recursive models to assist 
the reader in recognizing models which can be represented by 
the graphical structure. 

In 1 17 1 Shi and Nehorai consider a model for recursive scat- 
tering between point scatterers in a homogeneous background. 
The propagation between any point in space is described by 
a scalar Green's function. The transfer function obtained by 
applying the Foldy-Lax equation can also be obtained from 
a propagation graph by defining the sub-matrices of A(/) as 
follows. The model does not include a directed term and thus 
D(/) = 0. The entry of [T(/)]mj„ is the Green's function 
from transmit vertex mi to scatterer n times the scattering 
coefficient of scatterer n'. Similarly, the entry [R(/)]„,„2 is 
the Green's function from the position of scatterer n to receiver 
1712- The entry [B(/)]„„',n ^ n' is the Green's function from 
the position of scatterer n to the position of scatterer n' times 
the scattering coefficient of scatterer n. Since a point scatterer 
does not scatter back on itself, the diagonal entries of B(/) 
are all zero. As can be observed from the above definitions, 
the assumption of homogeneous background medium leads to 
the special case with 8d — %, £t — Vt^ Vs, = Vs x Vr, and 

= Vs X Vs. 

Another modeling method that can be conveniently de- 
scribed using propagation graphs is (time-dependent) radiosity 
lfT9l -ll23l. In these methods, the surfaces of the objects in the 
propagation environment are first divided in smaller patches, 
then the power transfer between pairs of patches is defined 
in terms of the so-called "delay-dependent form factor", and 
finally the power transfer from the transmitter to the receiver 



is estimated. The delay-dependent form factor is essentially 
a power impulse response between patches. It appears that 
for these algorithms no closed-form solution feasible for 
numerical evaluation is available in the hterature. Thus |fT9l - 
p3l resort to iterative solutions which can be achieved after 
discretizing the inter-patch propagation delays. The time- 
dependent radiosity problem can be expressed in the Fourier 
domain in terms of a propagation graph where each patch 
is represented by a scatterer, while the entries of A(/) are 
defined according to delay dependent form factor Using this 
formulation, a closed form expression of the channel transfer 
matrix appears immediately by Theorem [T] with no need for 
quantization of propagation delays. 

F. Revisiting Existing Stochastic Radio Channel Models 

It is interesting to revisit existing radio channel models by 
means of the just defined framework of propagation graphs. 
Such an effort may reveal some structural differences between 
models, which are not apparent merely from the mathematical 
formulation. It is, however, a fact that the interpretation of 
a transfer function as a propagation graph is not unique — 
it is straightforward to construct different propagation graphs 
which yield the same transfer matrix. Thus different propa- 
gation graphs may be associated to the same radio channel 
model. In the sequel we construct graphs for two well-known 
stochastic channel models, i.e., the model by Turin et al. Il33l 
and the model by Saleh and Valenzuela fS?]. We include here 
only enough detail to allow for the discussion of the structure 
of the associated graphs and refer the reader to the original 
papers for the full description. 

The seminal model [33,1 by Turin et al. can be expressed in 
the frequency domain as 

oo 

H{f ) = E exp(-j2^/Tf), (29) 

£=0 

where is the complex gain and ti denotes the delay of 
component £. Apart from the direct component, to which we 
assign index £ = 0, the index £ is merely a dummy index 
which does not indicate any ordering of terms in (|29^ . The set, 
{(rf, ag) : £ = 1, 2, . . . } is a marked Poisson point process on 
[tq, cxd) with complex marks {ag : I — 1,2, . . .}. 

The propagation graph we seek should contain a propaga- 
tion path for each term in the sum (l29T l. We notice that mod- 
ifying the term exp(— j'27r/T^) has no effect on the other 
terms in (|29] |. This structure can be captured by constructing 
a propagation graph in which modifying vertices or edges in a 
propagation path leaves others unchanged. Graph theory offers 
the convenient notion of "independent walks" ifJTI : Two or 
more walks are called independent if none of them contains 
an inner vertex of another In particular, propagation paths are 
independent if they do not traverse the same scatterers. Thus 
for independent propagation paths changing path £ by 
changing its edge transfer functions or by deleting edges from 
it, has no effect on path £' . We therefore associate to each 
term in (|29] | an independent propagation path. Assuming for 
simplicity all paths to be single bounce paths this results in 
the graph shown in Fig. |3a). 
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a) Turin et al. 



Tx ■ 




Rx 



the experimentally observed avalanche effect and diffuse tail 
can be explained using only scatterer interactions that by them- 
selves yield no delay dispersion. We specify a method feasible 
for generating such a graph in Monte Carlo simulations. The 
model discussed in this example is a variant of the model 
proposed in ||29l. Il30l. 



b) Saleh & Valenzuela 



Tx ■ 





Rx 



Fig. 5. Propagation graph representations of: a) a realization of the model 
by Turin et al. 1331 ; and b) a realization of the Saleh- Valenzuela model 1341 . 



The model by Saleh and Valenzuela Il34l can be formulated 
as a second-order Turin model: 

oo oo 

H{f ) = exp(-i2^/T,) aw exp(-j2^/T«0, (30) 



where ^ is a cluster index and £' is a index for the components 
within a cluster The set {r^ : £ ~ 1,2,...} is a Poisson 
process on [ro,oo) and conditioned on the cluster delays 
{Te;£ = 0,1,...}, the family of sets {{ctu',Tu') ■ = 
1,2,...}, £ — 0,1,..., are independent marked Poisson 
processes. Using a similar argument as for the Turin model, 
leads to the graph depicted in Fig. lUb). The transmitter is 
connected to a set of super-ordinate scatterers, each scatterer 
corresponding to a "cluster". These cluster-scatters are then 
connected to the receiver via independent single-bounce paths 
passing through sub-ordinate scatterers. In this graph, deleting 
cluster scatterer £ makes a whole signal cluster £ disappear 
in dSOl l. Similarly, removing the sub-ordinate scatterer £' the 
component with index ££' in (|30] | vanishes. 

We end the discussion of the Saleh- Valenzuela model by 
mentioning that many equivalent graphical interpretations may 
be given. As an example, one may consider a graph structure 
as that of the reverse graph of Fig. 0b). In such a structure, the 
transmitter is directly connected to the sub-ordinate scatterers 
whereas the receiver is connected to the clusters. Indeed by 
reversing any number of clusters we obtain an equivalent 
propagation graph. These graphs share the common property 
that they contain no cycles. 

IV. Example; Stochastic Model for In-Room 
Channel 

The concept of propagation graph introduced until now 
can be used for describing a broad range of channel models. 
In this section we apply these general results to a specific 
example scenario where scatterer interactions are considered 
to be non-dispersive in delay. This resembles the case where 
all wave interactions can be considered as specular reflections, 
or considering points scatterers. The intent is to exemplify that 



A. Weighted Adjacency Matrix 

We define the weighted adjacency matrix according to a 
geometric model of the environment. We consider a scenario 
with a single transmitter, a single receiver, and Ng scatterers, 
i.e., the vertex set reads V = VtUVi UVs with Vt = {Tx}, Vr = 
{Rx}, and Vs = {SI, . . . , SiVs}- To each vertex v e V 
we assign a displacement vector r„ £ M.^ with respect to a 
coordinate system with arbitrary origin. To edge e — {v, v') 
we associate the Euclidean distance = ||r„ — r^,/ 1|, the gain 
ge, the phase 0e, and the propagation delay Tg — d^/c, where 
c is the speed of light. The edge transfer functions are defined 
as 



AM) 



9e{f)expj{ 
0; 



27rre/); e € £ 



(31) 



The edge gains {ge{f)} are defined according to 



g . 

odi(e)^ ' 



e e f d 
e e £t 
e e Er 
e e £s 



(32) 



where odi(e) denotes the number of edges from init(e) to 
other scatterers and for any £' C £ 



and 



Si£') = 



eeS' 



ee£' 



(33) 



with I • I denoting cardinality. The weight of the dkect edge 
is selected according to the Friis equation ||36l assuming 
isotropic antennas at both ends. The weights of edges in £t 
and £i- also account for the antenna characteristics. They are 
computed at the average distance to avoid signal amplification 
when scatterers are close to a transmitter or receiver, namely 
when the far-field assumption is invalid. The value of inter- 
scatterer gain g is for simplicity assumed fixed. 

B. Stochastic Generation of Propagation Graphs 

We now define a stochastic model of the sets {r„}, £, and 
{0e} as well as a procedure to compute the corresponding 
transfer function and impulse response. The vertex positions 
are assumed to reside in a bounded region 7?. C K."^ corre- 
sponding to the region of interest. The transmitter and receiver 
positions are assumed to be fixed, while the positions of the 
Ns scatterers {r„ : f G Vs} are drawn independently from a 
uniform distribution on TZ. 

Edges are drawn independently such that a vertex pair e S 
V X V is in the edge set £ with probability Pg = Pr[e e £] 



g 



1) Draw r„, v £Vs uniformly on TZ 

2) Generate £ according to ( I34t 

3) Draw independent phases {^e '■&££} uniformly on 

[0,2^) 

4) Compute A(/) within the frequency bandwidth using J31t 

5) IF spectral radius of B(/) is larger than unity for some 
frequency within the bandwidth THEN GOTO step 1 

6) Estimate HK-.L^f) and hK:L{T) as described in AppendiXp 

Fig. 6. Algorithm for generating full or partial transfer functions and impulse 
responses for a preselected bandwidth. 

defined as 



>dir, 


e = (Tx, Rx) 


0, 


terni(e) = Tx 


< 0, 


init(e) = Rx 


0, 


init(e) — term(e) 


^ VIS 


otherwise 



The first case of ( l34l i controls the occurrence of a direct 
component. If Pdir is zero, D{f) is zero with probability 
one. If Pdir is unity, the direct term D{f) is non-zero with 
probability one. The second and third cases of (l34t exclude 
ingoing edges to the transmitter and outgoing edges from the 
receiver. Thus the generated graphs will have the structure 
defined in Section III-BI The fourth case of ( l34b excludes the 
occurrence of loops in the graphs. This is sensible as a specular 
scatterer cannot scatter a signal back to itself. A consequence 
of this choice is that any realization of the graph is loopless 
and therefore A(/) has zeros along its main diagonal. The 
last case of ( l34l i assigns a constant probability Pvis of the 
occurrence of edges from Vt to Vs, from to and from 

Vs to Vr. 

Finally, the phases {(/)e : e G f } are drawn independently 
from a uniform distribution on [0; 2tt). 

Given the deterministic values of parameters TZ, rxx, trx, 
-^s, ^'dii, ^vis and g, realizations of the (partial) transfer func- 
tion HK:L{f) and corresponding (partial) impulse response 
hK:L{T) can now be generated for a preselected frequency 
range [/min,/max], using the algorithm stated in Fig. |6] and 
Appendix. 

C. Numerical Experiments 

The effect of the recursive scattering phenomenon can 
now be illustrated by numerical experiments. The parameter 
settings given in Table|T]are selected to mimic the experimental 
setup of li26l used to acquire the measurements reported 
in Fig. [T] The room size and positions of the transmitter 
and receiver are chosen as in ll26l . We consider the case 
where direct propagation occurs and set Pdir to unity. The 
probability of visibility Pvis and the number of scatterers 
A^s are chosen to mimic the observed avalanche effect. The 
value of g is set to match the tail slope p ?» — 0.4dB/ns 
of the delay-power spectrum depicted in Fig. [T] The value 
of g can be related to the slope p of the log delay-power 
spectrum via the approximation p k, 20\ogif^{g) / p{£s). This 



TABLE I 

Parameter Settings for the Numerical Examples 



Parameters 


Symbol 


Values 


Room size 


TZ 


^U, oj X [U, oj X [u, 2.o\ m 


Transmitter position 




[1.78, 1.0, 1.5]'^m 


Receiver position 




[4.18, 4.0, 1.5]'^ m 


Number of scatterers 




10 


Tail slope 


P 


-0.4dB/ns 


Prob. of visibility 




0.8 


Prob. of direct propagation 


Pdir 


1 


Speed of light 


C 


3-10" m/s 


Transmit signal 


X[m] 


Unit power Hann pulse 


Number of frequency samples 


M 


8192 




2 2.2 2.4 2.(1 2.8 3 




TT^,n:^20 40 60 80 100 120 

T [lis] 



Fig. 7. Channel response for a specific realization of the propagation graph. 
Top: Transfer function in dB (201ogio \H{f)\). Bottom: Impulse responses 
in dB (201og]^Q |/i(r)|) computed for two frequency bandwidths. 



approximation arises by considering the power balance for a 
scatterer: assuming the signal components arriving at the scat- 
terer to be statistically independent, neglecting the probability 
of scatterers with outdegree zero, and approximating the delays 
of edges in Eg by the average p{£s) defined in ( l33T l. 

Fig. |7] shows the amplitude of a single realization of the 
transfer function. Overall, the squared amplitude of the transfer 
function decays as /^^ due to the definition of {gdf)}- 
Furthermore, the transfer function exhibits fast fading over the 
considered frequency band. The lower panel of Fig. |7]reports 
the corresponding impulse response for two different signal 
bandwidths. Both impulse responses exhibit an avalanche 
effect as well as a diffuse tail of which the power decays 
nearly exponentially with p « — 0.4dB/ns. As anticipated, 
the transition to the diffuse tail is most visible in the response 
obtained with the larger bandwidth. 

The build up of the impulse response can be examined via 
the partial impulse responses given in Fig. |8] Inspection of the 




Fig. 8. Partial responses obtained for one graph realization for the bandwidth [2, 3] GHz. The K : L settings are indicated in each miniature. Top row: 
-bounce approximations. Right-most column: en'or terms resulting from {K — l)-bounce approximations. Main diagonal {K = L): _R'-bounce contributions. 



partial responses when K = L reveals that the early part of 
the tail is due to signal components with a low K while the 
late part is dominated by higher-order signal components. It 
can also be noticed that as K increases, the delay at which 
the maximum of the X-bounce partial response occurs and 
the spread of this response are increasing. 

Fig. |9] shows two types of delay-power spectra. The upper 
panel shows the ensemble average of squared amplitudes of 
1000 independently drawn propagation graphs for the two 
signal bandwidths also considered in Fig. |7] Both spectra 
exhibit the same trend: A clearly visible peak due to the direct 
signal is followed by a tail with nearly exponential power 
decay. As expected, the first peak is wider for the case with 
1 GHz bandwidth than for the case with 10 GHz bandwidth. 
The tails differ by approximately 7 dB. This shift arises due 
to the /^^ trend of the transfer function resulting in a higher 
received power for the lower frequencies considered in the 
1 GHz bandwidth case. 

The bottom panel shows spatially averaged delay-power 
spectra obtained for one particular realization of the propaga- 
tion graph. The simulated spatial averaged delay-power spectra 
exhibit the avalanche effect similar to the one observed in 
Fig. [T] For the 10 GHz bandwidth case the power level of the 
diffuse tails agrees remarkably well with the measurement in 




[1.11] [GHz] 



[1.11] [GHz] 




Fig. 9. Simulated delay-power spectra. Top panel: Ensemble average of 1000 
Monte Carlo runs. Bottom panel: Spatial average of a single graph realization 
assuming the same grid as in Fig. [T] i.e., 900 receiver positions on a 30x30 



horizontal square grid with 1x1 cm 
in Table H] 



mesh centered at position r^x given 
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Fig. [T] The modest deviation of about 3 dB can be attributed 
to antenna losses in the measurement. 

V. Conclusions 

The outset for this work was the observation that the mea- 
sured impulse responses and delay-power spectra of in-room 
channels exhibit an avalanche effect: Multipath components 
appear at increasing rate and gradually merge into a diffuse tail 
with an exponential decay of power We considered the ques- 
tion whether the avalanche effect is due to recursive scattering. 
To this end we modelled the propagation environment as a 
graph in which vertices represent transmitters, receivers, and 
scatterers and edges represent propagation conditions between 
vertices. From this general structure the graph's full and partial 
frequency transfer matrices were derived in closed form. These 
expressions can, by specifying the edge transfer functions, be 
directly used to perform numerical simulations. 

We considered as an example a graph-based stochastic 
model where all interactions are non-dispersive in delay in 
a scenario similar to an experimentally investigated scenario 
where the avalanche effect was observed. The impulse re- 
sponses generated from the model also exhibit an avalanche 
effect. This numerical experiment lead to the observation 
that the diffuse tail can be generated even when scatterer 
interactions are non-dispersive in delay. Furthermore, the ex- 
ponential decay of the delay-power spectra can be explained 
by the presence of recursive scattering. As illustrated by the 
simulation results the proposed model, in contrast to existing 
models which treat dominant and diffuse components sepa- 
rately, provides a unified account by reproducing the avalanche 
effect. 

Appendix 

The transfer function H(/) and impulse response h(r) can 
be estimated as follows: 

1) Compute M samples of the transfer matrix within the 

bandwidth [fm\n, /max] 

HH =H(/„,i„ + mA/), m = 0,l,...M-l, (35) 

where A/ = (/max - /min)/(M - 1) and H( ■ ) is 
obtained using Theorem [T] 

2) Estimate the received signal y(T) via the inverse discrete 
Fourier transform 

M-l 

y(iA^) = A/ ^ H[TO]X[m] exp(j2 

m=0 

i = Q,...M -I, 
where X[m] = X(/min + mAf ),m = 0, 1, . . . Af - 1 

and A^ = l/(/max - /min)- 

The impulse response can be estimated by letting X[/] be a 
window function of unit power 

/ |X(/)|2d/« ^ IXHII'A/ = 1, (36) 

"'/min ,„ = 

where X[to] = X(/nii„ + niAf). The window function must 
be chosen such that its inverse Fourier transform exhibits a 



narrow main-lobe and sufficiently low side-lobes; y(r) is then 
regarded as a good approximation of the impulse response 
of the channel and by abuse of notation denoted by h(r). 
Samples of the partial transfer matrix are obtained by replacing 
H( • ) by Hx:l( ■ ) in (l35T l. The corresponding received partial 
impulse response is denoted by hif:i(r). 
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